Spinodal of supercooled polarizable water 
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We develop a series of molecular dynamics computer simulations of liquid water, performed with a 
polarizable potential model, to calculate the spinodal line and the curve of maximum density inside 
the metastable supercooled region. After analysing the structural properties, the liquid spinodal line 
is followed down to T = 210 K. A monotonic decrease is found in the explored region. The curve of 
maximum density bends on approaching the spinodal line. These results, in agreement with similar 
studies on non polarizable models of water, are consistent with the existence of a second critical 
point for water. 
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I. INTRODUCTION 

Liquid water behaviour as a function of temperature 
and pressure differs from that of most liquids. As a major 
example the presence of a line of maximum density (TMD 
line) in the P-T phase diagram is most significant, as it di- 
vides the entire P-T phase diagram into two regions with 
remarkably different properties: the coefficient of thermal 
expansion is negative on the low-temperature side of the 
TMD line, while it is positive on the high-temperature 
side. Similarly, several other thermodynamic and dy- 
namical quantities show peculiar properties. For instance 
the isothermal compressibility and heat capacity exhibit 
a minimum as a function of temperature and show an 
anomalous behavior more pronounced in the supercooled 
liquid state When decreasing temperature in fact 
the coefficient of thermal expansion, the isothermal com- 
pressibility and the constant-pressure specific heat in- 
crease rapidly. These quantities appear to diverge at a 
temperature T w —45 C if they are extrapolated below 
the lowest temperatures at which they are measurable 
(T = -42 C at P = 1 bar) The supercooled region 
below T = —42 C is experimentally unreachable due to 
the strong tendency of water to crystallize. It is nonethe- 
less believed that nucleation in supercooled water might 
be due to the presence of impurities that drive the liquid 
toward the more stable phase [2|]. Therefore the exper- 
imentally unreachable zone of supercooling is believed 
to be physically significant. Measurements of the rate of 
evaporation on amorphous water have in fact proved that 
the amorph can be connected with normal liquid water 
by a reversible thermodynamic path at atmospheric pres- 
sure Q. Besides at ambient pressure experiments have 
proven the existence of supercooled liquid water close to 
the glass transition temperature [3]. 
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In spite of all the interest driven by these anomalies 
a coherent theory of the thermodynamic and transport 
properties of supercooled water does not yet exist, also 
due to the difficulties encountered in experiments. 

Different thermodynamic scenarios have been proposed 
through the years for the peculiar metastable behaviour 
of water [j| . Among these three of them have received a 
great attention in literature. The stability limit conjec- 
ture (SLC) is the first scenario proposed. It attributes 
the anomalies of water to the presence of a continuous 
retracing spinodal curve, bounding the superheated and 
supercooled states. The spinodal, which represents the 
limit of mechanical stability of a liquid, is hypothized 
to retrace to higher pressure values below a temperature 
at which it intersects the locus of the TMD @, l|. The 
second critical point scenario (SCP), based on extrapola- 
tion of simulated data, ascribes the anomalous properties 
of water to the presence of a metastable, low tempera- 
ture liquid- liquid critical point, associated with a phase 
transition between a low-density and a high density liq- 
uid phase Q . Experimental evidences of the existence of 
two liquid phases compatible with a second critical point 
have been presented [9(. The singularity free scenario 
(SF) explains the behaviour of supercooled water with 
the presence of anomalous fluctuations due to the hydro- 
gen bonds and no underlying singularity is invoked [To| . 
The SF assumes that upon isobaric cooling the thermo- 
dynamic response functions go through a maximum but 
remain finite. An anomalous increase in isothermal com- 
pressibility, heat capacity, and thermal expansion is ex- 
plained by the existence of the density maxima locus, 
which is negatively sloped in the P, T plane [2j . Accord- 
ing to this scenario, no phase transition or critical point 
occurs at low temperatures. 

A fundamental role in the clarification of the scenario 
of supercooled water is played by the study of the be- 
haviour of the spinodal line, and its relation with the 
line of TMD. Due to the limitations of the experiments in 
the supercooled realm, in order to assess the different hy- 
pothesis numerical studies have become of uttermost im- 
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portance. Many computer simulation studies indeed have 
been performed with several different water site models, 
ST2,TIP4P,SPC/E and T1P5P 0, d M, 0, M, d- 
All the above models do not take into account explic- 
itly the polarizability of the water molecules. Polarizable 
models for water have been also developed [13, d, 
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realistic description of the behaviour of the system in 
different thermodynamic conditions. Depending on the 
zone of the phase diagram investigated they have proven 
to be equally or more realistic to describe the features of 
water. Studies of the behaviour of the spinodal line for 
polarizable water model potential have been never car- 
ried out. It is therefore of interest to study the behaviour 
of the spinodal line for potentials where polarizability is 
explicitly taken into account. 

We present here a Molecular Dynamics, MD, study of 
the spinodal and the TMD line of supercooled water per- 
formed with the polarizable potential, introduced in the 
literature by Ruocco and Sampoli [22j and parametrized 
by Brodholt, Sampoli and Vallauri [23], the BSV po- 
tential. We found that the BSV model is appropriate 
for our study since it reproduces the thermodynamical 
properties of water in the region of the gas-liquid co- 
existence better than other polarizable models over a 
broad temperature range (30j |. It has also been shown 
that within the BSV model the site-site radial distribu- 
tion functions are in good agreement with the experi- 
mental data in a broader range of thermodynamical con- 
ditions, [m, m m, m, m 

In the next section we describe the computer simula- 
tion in details. After presenting in the third section the 
structural properties of the system upon supercooling, 
the behaviour of the spinodal and of the TMD lines are 
reported in the fourth section. Last section is devoted to 
conclusions. 



II. COMPUTER SIMULATIONS 

We performed MD computer simulations of water with 
the Model 4 of ref. 23]. An induced polarizable dipole 
moment P, located in the center-of-mass of the molecule, 
describes the effect of the electric field of the environment 
on the molecules. The induced dipole pi — aEi is cal- 
culated from the local electric field Ei with an iterative 
procedure by assuming an isotropic polarizability. The 
value of this polarizability is fixed to the single molecule 
value a = 1.44 A 3 

Each simulation is conducted in the NVT ensemble 
with 256 water molecules enclosed in a cubic box with pe- 
riodic boundary conditions. The simulations have been 
performed with the minimum image convention and a 
cut-off of the interactions at half of the box length. The 
long range part of the electrostatic interactions is taken 
into account with the reaction field method. Details of 
the extension of the reaction field method to include po- 
larization effects can be found in references [22, [39] . The 



time step 5t for the integration of the molecular trajec- 
tories is fixed at 1 fs. 

To analyse the phase diagram and the spinodal curve 
inside the metastable supercooled region we carried out 
isothermal simulations for eight different temperatures: 
350, 300, 280, 260, 240, 230, 220 and 210 K. For each 
eight different densities are simulated, namely 1.05, 1.00, 
0.98, 0.95, 0.90, 0.87, 0.85 and 0.83 gr/cm 3 . The box 
length L spans from L = 19.39 A to L = 20.96 A to 
cover the range of investigated densities. 

During equilibration the temperature is controlled via 
a velocity rescaling procedure. Production runs are per- 
formed in the microcanonical ensemble. For the temper- 
atures and the densities closer to the spinodal line the 
longest equilibration time required was of 5 ns. 

In t ables HI and [TT1 we report the temperatures, densities, 
pressures and internal energies of the simulations we have 
performed. 



III. STRUCTURAL PROPERTIES 

Structural quantities have been calculated averaging 
from a minimum of 2.5 x 10 3 configurations at T = 300 K 
to a maximum of 3 x 10 3 configurations at T = 220 K, all 
equally spaced and taken every 2ps. We have calculated 
the site-site radial distribution functions (RDF) g a p{r) 
to monitor the internal structure at every temperature. 
We report goo{ r ): gouir) and guH{r) for T = 300 K 
in FigEDand goo{r) for T = 220 K in FigH calculated 
at several densities. At T = 300 K the peak positions 
agree with the experimentally measured structure of wa- 
ter 37J. Upon decreasing temperature the peaks remain 
at the same position and become sharper. By varying 




FIG. 1: Site-site pair correlation functions of the BSV water 
at T = 300 K and densities p = 1.00, 0.98, 0.95, 0.90, 0.87, 
0.85, 0.83. goo(r) is in the main frame. The insets show 
goH{r) and gHn(r) starting from the top. 

the densities we do not observe very large changes in the 
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9a.p{r) apart for an increase in the height of the peaks, 
even when the system approaches the spinodal line (see 
next section) . This is an important test to verify the ab- 
sence of disomogeneities or cavitation phenomena close 
to the limit of mechanical stability. A further check has 
been performed with the calculations of the local density 
fluctuations. For this sake we have divided our simula- 
tion box into 64 sub-boxes and calculated the average 
distribution functions of the local density (not shown). 
In all cases we do not find deviations from a gaussian 
distribution centered on the fixed density of the system. 

A deeper insight in the local ordering of the atoms can 
be achieved by calculating the coordination numbers. In 
the upper inset of Fig. [2] we show noo , the coordination 
number of the oxygens, as a function of p: 



noo = 47T/7 



goo{r)r 2 dr 



(1) 



where r m i n is the value of the interatomic distance at 
which the first minimum in goo( r ) is located. It is im- 
portant to observe the different behaviour of noo at the 
two temperatures T = 300 K and T = 220 K on lowering 
the density. At T = 300 K for the highest density noo 
results to be 5.2. This value is the signature of an open 
network of water, where also interstitial molecules are 
present [33| . noo then decreases until it reaches the lim- 
iting value 4 characteristic of the local tetrahedral order. 
At the lower temperature T — 220 K noo is found to be 
quite constant around the value 4. These findings indi- 
cate that the open network of f/-bonds becomes better 
defined as p and/or T decrease. 

In connection with the tendency of the _ff-bond net- 
work to become more ordered it is expected a sharpen- 
ing of peaks of the RDF. A quantitative measure of this 
sharpening can be obtained looking at the value of the 
RDF at its first minimum, g(r m i n ). As a consequence we 
also expect to observe a decrease of g[r m in) as the coor- 
dination shells become better defined. In the lower inset 
of Fig. [2] we observe instead the presence of a minimum 
in the behaviour of g(r m in) as function of the density. 
The minimum is more evident at T = 220 K. This effect 
evidences an important change in the thermodynamical 
properties of the system, since it is related to the be- 
haviour of the entropy in the vicinity of the spinodal, as 
previously discussed in literature [1 1] . Therefore we can 
consider that the change of slope in the g(r m i n ) curve is 
a signature of the approaching of the system to the limit 
of mechanical stability. 



IV. SPINODAL LINE 

The limit of the mechanical stability of a system can 
be determined from a study of the isothermal compress- 
ibility. In the space of variables identifying the thermo- 
dynamic state of the system in fact, the condition for the 
mechanical stability is [|J 
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FIG. 2: Pair correlation function goo(r) at T = 220 K and 
the same densities as Fig. [T] In the inset on the top number 
of nearest neighbours noo as a function of p for T = 300 K 
(full squares) and T = 220 K (full circles); in the inset below 
value of the first minimum of goo{r), g{r m in), as a function 
of p at T = 300 K (full squares) and T = 220 K (full circles). 



where Kt is the isothermal compressibility defined as 



K T = - 



1 f dp 



p \dPJ T 



(3) 



K T >0 



(2) 



In the phenomenological description of the liquid-gas 
transition, as given for instance by mean field theories like 
the Van der Waals equation, along an isotherm the limits 
of the mechanical stability are marked by the changes in 
the slope of the Pt(p) curve. Therefore from Eq. {3J it is 
found that the crossing to unstable states with Kt < 
takes place at singularity points where Kt diverges. The 
spinodal line is identified by these singularity points. In 
particular the point where the minimum of Pt{p) is lo- 
cated represents the limit of stability of the superheated 
liquid and so it belongs to the liquid branch of the spin- 
odal. 

We are interested here at the liquid spinodal which 
extends in the P — T plane below the liquid-gas coexis- 
tence into the region of negative pressures. The liquid 
spinodal indicates also the limiting values of tension for 
the existence of an homogeneous fluid before phenomena 
like cavitation take place [|J • 

The behaviour of Pt{p) for the BSV water is shown 
in Fig. [3]for p ranging from 1.05 to 0.83 gr/cm 3 . In the 
inset of Fig. [3] four isotherms are represented upon low- 
ering the temperature from T — 350 K to T = 260 K. 
Isotherms from T = 240 K to T = 210 K are shown in 
the main frame of Fig. [31 From T — 300 K to the lowest 
temperature investigated the Pt(p) curves exhibit a min- 
imum in the simulated range of densities. The isotherms 
corresponding to lower temperatures have larger fluctu- 
ations due to the influence of both the low temperatures 
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and the negative pressures In order to better follow the 
curves for the lowest temperatures shown in Fig. [3] we fit- 
ted the calculated points of the isotherms with a fourth- 
order polynomial function. We observe however that due 
to the fluctuations the fit systematically slightly over- 
estimates the values of the minima. Therefore we will 
consider the minima directly observed for the plot of the 
spinodal line. 

Based on Eq. [2] the minima of the isotherms represent 
the limit of stability of the superheated liquid and iden- 
tify the values of the pressure at the spinodal, P S {T). 
As T decreases P S (T) clearly shifts to larger values o f p . 
This behaviour is similar to that observed for TIP5 P lla ] 
while for ST2 and TIP4P this shift is less evident [ljj7^ 

The resulting estimated liquid spinodal line P S (T) is 
plotted in fig. 2] In the same figure also the P P (T) iso- 
chores are shown as obtained by the data set of tables U 
andflj] At least down to the lowest investigated tempera- 
ture, T = 210 K, the spinodal line is not retracing. In the 
SLC framework the spinodal line is expected to become 
reentrant at the point where it intersects the TMD line 
in the phase diagram. We recall that this behaviour is 
required if the TMD line has a negative slope @, [HI [3!| . 

Along the TMD line the coefficient of thermal expan- 
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FIG. 3: Isotherms (symbols) for densities ranging from 1.05 
to 0.83 gr/cm 3 . Errors are calculated with the method of 
statistical inefficiency [38J. In the inset temperatures range 
from 350 K to 260 K. In the main frame temperatures range 
from 240 K to 210 K and lines are fit with a fourth order 
polynomial function. 
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goes to 0. Since the thermal pressure coefficient 
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df 



(4) 



(5) 



is connected to the coefficient ap by the following relation 



Iv 



ap 



(6) 



the TMD points are on the line connecting the minima 
of the estimated P P {T) isochore for different p. 

According to Fig. [4] as p decreases, approaching the 
spinodal line, the isochores exhibit a minimum at smaller 
values of T. 

The isochore for the lowest reported density p = 
0.90 gr/cm 3 does not show any minimum in the range 
of T investigated where instead the minima for the other 
isochores are found. The TMD line does not approach 
the spinodal line and at the lower T investigated it results 
to have a positive slope. This behaviour of the TMD line 
prevents an intersection with the spinodal at negative P 
and low T, unless a further change in the slope of the 
TMD line is hypothized. 

These results are in agreement with the findings of the 
studies done on non polarizable models for water [TJ [11 
El QUE OS- In fact for all these models the TMD line 
changes slope to avoid intersection with a non reentrant 
spinodal. Differences among the models can be found 
in the values of pressures and temperatures where the 
curves span. In particular the change of slope of the 



TMD line appears in our model at negative pressures, 
P < -120 MPa, similar to SPC/E [H|, P < -80 MPa 
and to ST2 (Tj. For TIP5P [Q the change of slope 
appears at ambient pressure instead. 

The most studied alternative to the SLC for explaining 
the anomalies of water is the presence of a second critical 
point. The SCP hypothesis H [l4| is based on the con- 
sideration that in the glassy state two distinct types of 
structure are observed [40], the low density amorphous 
(LDA) ice and the high density amorphous (HDA) ice. 
This polymorphism of water would be extended into the 
liquid phase and it would emerge as a coexistence line 
between two liquid phases culminating in a second crit- 
ical point [4l[. Computer simulation of waterlike lattice 
models show evidence in fact of two immiscible liquid 
forms of water El El Si El . 



More recently several amorphous ice phases [46| have 
been experimentally observed and signatures have been 
found in computer simulation of the existence of sev- 
eral liquid- liquid coexistences in supercooled water [47| . 
Also Gibbs ensemble Monte Carlo (GEMC) studies on 
the BSV water have found two liquid-liquid coexistence 
curves besides the liquid- vapor equilibrium (30j . The liq- 
uid branch of the gas-liquid coexistence curve obtained 
with the GEMC is reported in Fig. [5j where we collect 
the relevant features of the behaviour of the BSV water 
in the (T, p) thermodynamical plane. The spinodal line 
obtained in the present work in the supercooled region 
is compatible with the GEMC results. It comes out that 
the metastability region explored here goes very close to 
the boundaries of the region where new phases of the 
liquid water have been found. 
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FIG. 4: P P (T) isochores for several values of p in g/cm 3 , 
spinodal line and TMD curve. As p decreases, approaching 
the spinodal line, the isochores exhibit a minimum at smaller 
values of T. The spinodal line decreases upon lowering tem- 
perature and no retracing behaviour is observed. 
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FIG. 5: Liquid branch of the coexistence curve obtained by 
GEMC 30] (filled triangles) and the spinodal line of the 
present work (black circles) in the T, p plane for the BSV 
water. The diamond represents the estimated gas-liquid crit- 
ical point C. Portions of the liquid-liquid coexistence regions 
are shown by the open triangles [30l |. 



V. CONCLUSIONS 

We have done Molecular Dynamics simulation of su- 
percooled water described with a polarizable BSV poten- 



tial. The aim of the paper was to calculate the behaviour 
of the spinodal line and of the TMD line down to lowest 
possible temperatures and pressures. 

We have performed a structural analysis of water de- 
scribed by the BSV polarizable potential in the region of 
the supercooled liquid and we have found that the system 
persists in an homogeneous phase down to T = 210 K, 
P k, —390 MPa and this occurs for p m 0.9. Signatures 
of the presence of the mechanical stability limit are ob- 
served in the behaviour of the values of the RDF at the 
first minimum. 

The boundaries of the liquid mechanical stability at 
low T and negative P have been determined by the direct 
calculations of the liquid spinodal trough the minima of 
the isotherms Pr(p)- The line of the TMD in the same 
thermo dynamical region has also been detected from the 
minima of the isochores P p (T). 

We did not find any retracing of the spinodal line in 
the investigated region while we observe a change of slope 
of the TMD line which prevents for continuity a possi- 
ble intersection of the TMD with the spinodal even down 
to lower T and P. These results are in agreement with 
previous calculations done with non polarizable site po- 
tentials [H El El EI El El- Our analysis excludes 
indeed the possibility of the SLC scenario. 

As a matter of fact the no retracing of the spin- 
odal is compatible however with both the SI interpre- 
tation Jl(| and the possible existence of a second critical 
point On the other hand we have shown that the 
liquid spinodal calculated in the present work is in agree- 
ment with t he g as-liquid coexistence curve obtained with 
the GEMC Since in the GEMC there are clear in- 
dications of the existence of two liquid-liquid coexistence 
regions we can infer that our results support the second 
critical point scenario. 

The use of a polarizable potential in the present work 
has given an answer to two important issues. The first 
is that polarizable and non polarizable potentials are 
equally able to describe an unique, and therefore likely, 
scenario for the liquid-gas spinodal of supercooled water. 
This is an important achievement as the parameters of 
these potentials have been calculated through fit to ex- 
perimental data performed at ambient conditions. There- 
fore the uniqueness of phase behaviour far from these 
conditions reinforces the possibility that these potentials 
are able to predict existing features of water. The second 
is that the analysis performed with the polarizable model 
can be extended in a different zone of the phase diagram 
for determining the liquid-liquid spinodal and indications 
about the second critical point location. 
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350 


1.00 


100 ±2 


-39.16 ±0.17 


350 


0.98 


48 ± 24 


-39.03 ±0.18 


350 


0.95 


-23 ±4 


-37.63 ± 0.04 


350 


0.90* 


-112 ±40 


-37.28 ± 0.04 


350 


0.87 


-150 ±4 


-36.34 ±0.18 


350 


0.86* 


-159 ±4 


-36.65 ± 0.06 


350 


0.85 


-164 ±38 


-36.52 ± 0.06 


350 


0.83 


-182 ±4 


-35.58 ± 0.04 


300 


1.05 


121 ±5 


-42.92 ±0.13 


300 


1.00 


-22 ±3 


-41.55 ±0.10 


300 


0.98 


-72 ± 10 


-42.11 ± 0.05 


300 


0.95 


-128 ±20 


-40.97 ±0.02 


300 


0.90* 


-207 ±37 


-40.76 ±0.15 


300 


0.87 


-240 ± 37 


-39.72 ± 0.04 


300 


0.86* 


-243 ± 25 


-40.02 ± 0.07 


300 


0.85 


-238 ± 5 


-39.49 ± 0.04 


300 


0.83 


-230 ± 6 


-39.02 ± 0.16 


280 


1.05 


85 ±4 


-44.10 ±0.14 


280 


1.00 


-78 ±4 


-43.25 ± 0.04 


280 


0.98 


-112 ±4 


-43.77 ±0.08 


280 


0.95 


-170 ±6 


-42.56 ± 0.06 


280 


0.90* 


-242 ± 34 


-42.77 ±0.06 


280 


0.87 


-274 ± 5 


-41.62 ±0.09 


280 


0.86* 


-280 ± 4 


-41.71 ±0.10 


280 


0.85 


-281 ± 30 


-41.06 ± 0.12 


280 


0.83 


-234 ± 26 


-41.57 ±0.11 


260 


1.05 


42 ±36 


-45.90 ± 0.10 


260 


1.00 


-97 ±7 


-44.91 ± 0.03 


260 


0.98 


-132 ±5 


-45.67 ±0.11 


260 


0.95 


-200 ± 7 


-44.47 ± 0.08 


260 


0.90* 


-285 ± 7 


-44.65 ±0.10 


260 


0.87 


-312 ±27 


-44.26 ± 0.08 


260 


0.86* 


-300 ±21 


-44.18 ±0.09 


260 


0.85 


-270 ± 30 


-44.38 ± 0.06 


260 


0.83 


-216 ±28 


-43.90 ± 0.08 



TABLE I: Temperature, density, pressure and internal energy 
of the simulated state points. 
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T{K) 


P (9/ cm 3 ) 


P±SP (MPa) 


U±SU (KJ/mol) 


240 


1.05 


22 ± 10 


-47.09 ± 0.08 


240 


1.00 


-106 ±8 


-46.81 ± 0.06 


240 


0.98 


-163 ±6 


-47.48 ± 0.09 


240 


0.95 


-209 ± 28 


-46.52 ± 0.06 


240 


0.90* 


-305 ± 18 


-46.87 ±0.14 


240 


0.87 


-342 ± 25 


-45.82 ± 0.08 


240 


0.86* 


-329 ± 18 


-46.23 ± 0.06 


240 


0.85 


-268 ± 22 


-45.78 ± 0.07 


240 


0.83 


-236 ± 20 


-46.42 ± 0.08 


230 


1.05 


6±4 


-47.95 ± 0.03 


230 


1.00 


-135 ±5 


-48.26 ± 0.03 


230 


0.98 


-204 ± 27 


-47.87 ±0.09 


230 


0.95 


-217 ±8 


-47.62 ± 0.09 


230 


0.90* 


-325 ± 24 


-47.26 ± 0.09 


230 


0.87 


-336 ± 20 


-46.67 ±0.07 


230 


0.86* 


-330 ± 10 


-47.52 ± 0.10 


230 


0.85 


-337 ± 17 


-46.63 ± 0.08 


230 


0.83 


-257 ± 23 


-46.14 ±0.09 


220 


1.05 


-39 ± 33 


-49.67 ±0.08 


220 


1.00 


-90 ± 15 


-48.74 ±0.08 


220 


0.98 


-170 ±25 


-48.87 ±0.09 


220 


0.95 


-219 ±5 


-48.18 ±0.05 


220 


0.90* 


-335 ± 24 


-48.29 ±0.08 


220 


0.87 


-327 ±26 


-47.95 ± 0.05 


220 


0.86* 


-339 ± 14 


-47.78 ± 0.08 


220 


0.85 


-305 ± 11 


-46.85 ± 0.07 


220 


0.83 


-229 ± 26 


-47.38 ± 0.05 


210 


1.05 


9 ±30 


-50.09 ± 0.09 


210 


1.00 


-107 ±30 


-49.71 ±0.07 


210 


0.98 


-176 ±22 


-49.53 ± 0.08 


210 


0.95 


-226 ± 27 


-49.150 ±0.08 


210 


0.90* 


-393 ± 25 


-48.88 ±0.15 


210 


0.87 


-352 ±21 


-48.39 ± 0.09 


210 


0.86* 


-253 ± 30 


-49.10 ±0.11 


210 


0.85 


-263 ± 4 


-47.56 ± 0.08 


210 


0.83 


-337 ±24 


-47.86 ± 0.08 



TABLE II: Temperature, density, pressure and internal en- 
ergy of the simulated state points. 



